Distribution of Entropy Production in a Single-Electron Box 
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Recently, the fundamental laws of thermody- 
namics have been reconsidered for small systems. 
The discovery of the fluctuation re lations [3-0] 
has spurre d th eoretical [l|, @, @, l9 Ml3l | and experi- 
mental |l4t- |23| studies on thermodynamics of sys- 
tems with few degrees of freedom. The concept 
of entropy production has been extended to the 
microscopic level by considering stochastic trajec- 
tories of a system coupled to a heat bath. How- 
ever, the experimental observation of the micro- 
scopic entropy production remains elusive. We 
measure distributions of the microscopic entropy 
production in a single-electron box consisting of 
two islands with a tunnel junction. The islands 
are coupled to separate heat baths at different 
temperatures, maintaining a steady thermal non- 
equilibrium. As Jarzynski equality between work 
and free energy is not applicable in this case, the 
entropy production becomes the relevant param- 
eter. We verify experimentally that the integral 
and detailed fluctuation relations are satisfied. 
Furt herm ore, the coarse-grained entropy produc- 
tion |lol- 12. 23, 24] from trajectories of electronic 
transitions is related to the bare entropy produc- 
tion by a universal formula. Our results reveal 
the fundamental roles of irreversible entropy pro- 
duction in non-equilibrium small systems. 

Entropy production is a hallmark of irreversible ther- 
modynamic processes. The concept of a stochastic micro- 
scopic trajectory allows one to define entropy for small 
systems [if. However, such trajectories depend on the 
scale of observation. If one only accesses mesoscopic de- 
grees of freedom, one observes coarse-grained trajectories 
of mesoscopic states. The corresponding entropy produc- 
tion then differs from the bare entropy production with- 
out coarse-graining. In fact, it has been recently shown 
that coarse-graining of the slow background degrees of 
freedom for stochastic dynamics may actually lead to a 
modification of the fluctuation relations for entropy [13] ■ 
To clarify the concept of microscopic entropy production 



in non-equilibrium, accurate measurements are needed 
for systems, where the concepts of stochastic dynamics 
and time scale separation between the system and the 
heat bath are well-defined. 

A single-electron box (SEB) device at low tempera- 
tures is an excellent test bench for thermodynamics in 
small systems [12, [H, . 

The SEB employed here is shown in Fig. HJa). The 
electrons in the normal-metal copper island (N) can tun- 
nel to the superconducting Al island (S) through the alu 
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minum oxide insulator (I). The sample fabrication 
methods are similar to those in Ref. [22j, but the design 
is different in that the S side of the junction does not 
overlap with the normal conductor in order to intention- 
ally weaken the relaxation of energy in S [28| . Moreover, 
the main results in Ref. 221 were extracted from mea- 



surements at the temperature of 220 mK, whereas these 
measurements are conducted at 140 mK. Lower temper- 
ature further weakens the relaxation significantly [281 ] . 
leading to and elevated temperature in S. We denote by 
n the integer net number of electrons tunneled from S to 
N relative to charge neutrality. As we can monitor the 
charge state n with a nearby single-electron transistor 
(SET) shown in Fig. HJa), we take our classical system 
degree of freedom to be n. 

The device in Fig. UJa) can be represented with a clas- 
sical electric circuit, in which the energy stored in the 
capacitors and the voltage sources is given by [25 , 29|, 30[ 



H(n g ,n) = E c (n - n g f - e 2 n 2 J(2C g ), 



(1) 



where Eq is the characteristic charging energy, C g is the 
gate capacitance, n g = C g V g /e is the gate charge in units 
of the elementary charge e, and V g is the gate voltage 
which drives the system externally. Equation ((T|) gives 
the internal energy of the system. In an instantaneous 
single-electron tunneling event from n = k to n = k + 1 , 
the drive parameters stay constant and hence the work 
done to the system vanishes. Thus the first law of ther- 
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FIG. 1: (a) Sketch of the measured system together with a 
scanning electron micrograph of a typical sample. The colors 
on the micrograph indicate the correspondingly colored circuit 
elements in the sketch. The single electron box consists of 
a normal and a superconducting island that are connected 
through a tunnel junction. The charge state of the box is 
continuously probed with a SET detector, (b) Example trace 
of the measurement data under sinusoidal protocol for the 
drive V g , plotted in green. This trace covers three realizations 
of the forward process (V g from -0.1 to 1 mV), and three 
realizations of the backward process (V g from 1 to -0.1 mV). 
The SET current Idet, plotted in black, indicates the charge 
state of the box. The output of the threshold detection is 
shown in solid blue, with the threshold level as indicated by 
the dashed red line. The work and entropy production is 
evaluated for each realization separately, determined by the 
charge state trajectory. 



modynamics states that the generated heat is given by 

Q k = H(n g , fc) - H(n g , k + 1) = E c [2(n g - k) - 1]. (2) 

It has been recently demonstrated that when the SEB 
is in thermal equilibrium and the transition rates obey 
the detailed balance condition, the Jarzynski Equality 
(JE) (e~^ w ) = e~ /3AF relating the work done in the 
system W to its free energy change AF can be verified 
both theoretically and experimentally to a high degree 
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of accuracy 



30( | . In the present work, how- 
ever, the two environments consisting of the excitations 
in the normal metal and the superconductor are at differ- 
ent temperatures Tn = l/(fcs/3iv) and Ts = l/(fcs/?s), 
respectively, and hence the JE cannot be applied. Nev- 
ertheless, we expect that our system should obey the so- 
called intergal fluctuation theorem [l[ 



for the total entropy production As to t = As + As m 
given in terms of the increase of the system entropy 
As = \n{P[n(tf)}/P[n(0)]} and the medium entropy pro- 
duction As m . Here, P[n(t)] is the directly measurable 
probability of the system to be in state n at time in- 
stant t given the initial condition and the drive n g . We 
can express the entropy production of the medium as 
As m = PnQn + ftsQs, where Qn and Qs are the heat 
dissipated along the trajectory in the normal metal and 
in the superconductor, respectively. We can measure the 
total dissipated heat Q = Qn + Qs directly by monitor- 
ing n(t) with the SET and using Eq. ([2]). The only essen- 
tial assumption here is that the tunneling is elastic since 
the parameters of the Hamiltonian ([T]) can be measured 
independently. We can further obtain the conditional 
probability of Qn on Q b y so me additional assumptions 
(for technical details, see [27[) and hence the probability 
distribution of As m . 

On the other hand, medium entropy production can 
be defined by [3, [13 



A> 



r n _-m + (tj) 



(4) 



where the system is taken to make transitions at time 



instants t; 



to the state 



and 



■ % _-tn + (tj) and T n+ _t n _(tj) are the corresponding for- 



from the state n_ 

r n _-Mi+ {tj) and r„ 

ward and backward transition rates. We refer to As^ as 
the coarse-grained (cc) medium entropy production since 
it can be shown analytically (27j that in a single tunneling 
event 



-As 



-As 



;(Q) 



(5) 



where the average is taken over Qn for a fixed Q. This 
equality further implies (As m }g > As^(Q) and provides 
a physical interpretation for As^, the definition of which 
in Eq. Q coincides with the definition of medium entropy 
for general stochastic systems 13|. Note that by intro- 
ducing transition rates, we have implicitly assumed that 
the system is Markovian, a fact that can be experimen- 
tally verified in our setup. 

As mentioned above, we can obtain experimentally the 
probability distributions P f ±(As to t) and P^±(As^ t ) of 
Astot = As + As m and As^ t = As + As^, respec- 
tively, and hence access the integral fluctuation theorem 
of Eq. that should be satisfied by all the distribu- 
tions. Here, P_j. is the distribution for a forward driving 
protocol n Sj _j.(t) and P<_ corresponds to the backward 
protocol n g .^(t) = n g ,^(tf — t). In addition, we expect 
our system to satisfy so-called detailed fluctuation rela- 
tions 



P^(As t ot)/Pp±(-As tot ) 



(G) 



= e Ast ° 



-As., 



= 1. 



(3) 



In our experiments, we drive the system with the gate 
charge n g (t) = uq — Acos(irft), where no « A « 0.5. 
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Figure QJb) shows the applied drive and an example 
trace of the detector current. Clearly, two discrete cur- 
rent levels corresponding to the charge states n = and 
n = 1 are observable. Due to the low bath tempera- 
tures, 130 — 160 mK, the relatively high charging energy 
Ec ~ 162 fieV = 1.88 K x ks, and low driving frequen- 
cies / < 120 Hz, the system essentially always finds the 
minimum-energy state at the extrema of the drive. Thus 
we can partition the continuous measurement into legs 
of forward and backward protocols, for which the charge 
state and gate charge change from to 1 and 1 to 0, re- 
spectively. Conversion of the current trace from such a 
leg using threshold detection yields a realization for a sys- 
tem trajectory n(t) which is used in the ensemble average 
with unit weight to obtain the desired distributions. The 
charging energy, the temperatures of the normal metal 
and the superconductor, the tunneling resistance of the 
junction Rt ~ 1.7 Mfi, and the excitation gap of the 
superconductor A m 224 /ieV are determined experimen- 
tally [27 1 . 

Since the charge state corresponds to the ground state 
in the beginning and at the end of the drive, the system 
entropy change As in Eq. ([3]) and the free-energy change 
vanish, and we thus only need to obtain As m in order to 
assess if the fluctuation relations are satisfied. To deter- 
mine As™ in Eq. Q , the time-dependent tunneling rates 
rW 3 -(i) need to be measured. To extract them from an 
ensemble of forward and backward pumping trajectories, 
we show in Fig. [2ja) the observed probabilities for the 
system charge state to be n = 1 for the forward and back- 
ward drives denoted by Pl>. The rates can be solved by 
comparing the measured data to the outcome from the 
master equation [27J ■ An example of results obtained this 
way are shown in Fig. EJb). The rates from the standard 
sequential tunneling model [27j are in agreement with the 
experimentally obtained data. Here Tjy is assumed to be 
the temperature of the cryostat, while Tg is obtained for 
each measurement from the fit as listed in Table 1. 

The medium entropy production for a tunneling event 
with coarse graining, As"(n 9 ), extracted from the fitted 
rates is shown in Fig. [5] (c), demonstrating the signifi- 
cant effect of the overheating of the superconductor. For 
n g ea 0.4 — 0.6, the tunneling probability is primarily de- 
termined by the thermal excitations of the superconduc- 
tor and not by n g . Thus, the rates for different directions 
are almost equal and As™ is nearly vanishing. Only tun- 
neling events that occur outside this n g range contribute 
significantly to the cumulative entropy production. 

Table U presents a collection of exponential averages 
for entropy production as to test the integral fluctua- 
tion relation. Figure [3] shows the experimentally ob- 
tained distributions for work and entropy production 
with coarse graining together with the theoretical predic- 
tions for / = 40 Hz and different TV. For comparison, 
the prediction of JE at the bath temperature is shown by 
the dashed lines. As expected, the measured work distri- 



TABLE I: Measurement parameters and obtained averages 
for work and entropy production. Ts j0 /i is the S temperature 
matching the state n = 0/1, and W = W — AF 



Meas. 


/ 
(Hz) 


T N 
(mK) 


"0 


Ts,o 
(mK) 


2s, 1 
mK) 




) ( e - As - t ) 


(e _Astot ) 


1 


20 


130 


0.526 


174 


177 


93 


1.085 


1.063 


2 


40 


130 


0.516 


174 


177 


129 


1.064 


1.053 


3 


80 


130 


0.507 


176 


178 


180 


1.074 


1.083 


4 


20 


142 


0.513 


179 


181 


20 


1.064 


1.030 


5 


40 


142 


0.509 


179 


181 


30 


1.054 


1.047 


6 


80 


142 


0.505 


180 


181 


45 


1.096 


1.100 


7 


120 


141 


0.504 


181 


182 


68 


1.241 


1.324 


8 


40 


153 


0.502 


184 


184 


11 


1.095 


1.058 


9 


80 


153 


0.503 


184 


185 


15 


1.140 


1.139 


10 


120 


153 


0.502 


185 


186 


20 


1.301 


1.370 



butions in Fig. [3ja) do not follow from JE, which would 
assume just one temperature. The difference between Tn 
and Tg decreases with increasing Tn, and hence the dif- 




FIG. 2: (a) Probability for n — 1 as a function of time, ex- 
tracted from the measured ensemble of processes with drive 
frequency 40 Hz. The triangles pointing right indicate the 
probability for a forward process, and those pointing left in- 
dicate the probability for a backward process. The sinusoidal 
drive n g is displayed as a solid green line, (b) The tunneling 
rates obtained from measured probabilities [27| (shown in tri- 
angles) for measurements 1 — 3 and 8—10 listed in Table U 
with their corresponding fits (solid lines) . (c) As^ of Eq. Q 
for the transition — >• 1 as a function of n g , obtained from the 
tunneling rates in (b). The dashed lines show /3jvQ, which is 
equivalent to both As^; and As m if Ts = Tn- 
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FIG. 3: (a) P->{f3 N (W - AF)) distributions for 40 Hz for- 
ward process at different bath temperatures. The symbols 
show measured values (all panels), solid lines are numerical 
expectations (all panels), and dashed lines demonstrate what 
the distribution would be for Ts = Tjv. (b) Corresponding 
P_»(Astot) distributions, (c) P_>(Astot) distributions for sin- 
gle jump trajectories. 

ference between data and the dashed lines decreases as 
well. All the entropy distributions in Fig. Et^b), obtained 
from the same trajectories as the work distributions, sat- 
isfy the integral fluctuation theorem within the errors. 
The peaks in the distributions in the vicinity of zero arise 
from the relatively long time spent during the drive in the 
region n g ps 0.4 — 0.6, where the entropy production is 
nearly vanishing, see Fig. Hfc). 

Figure |4] (a) displays the probability distributions of 
entropy production at fixed TV = 142 mK for various 
frequencies. The tails of the distribution broaden and 
the peak at As££ t = sharpens with increasing frequency. 
Figurc@](b) shows the distributions for forward and back- 
ward processes. The distributions are overlapping, apart 
from the positions of the peaks near vanishing As^ t . The 
offset of the peaks is explained by different superconduc- 
tor temperatures for different tunneling directions lead- 
ing to ro->i(n g = 0.5) < ri_>o("-g = 0.5), and hence an 
offset of the point As^ t = away from n g — 0.5. Because 
the charge degeneracy point is the most probable point 
for the tunneling to take place, the peak is located there. 
For the n : — > 1 event, this corresponds to negative 
entropy production and positive production is observed 
for n : 1 — > 0. Different temperatures for different tun- 
neling directions can be justified by the difference in the 
observed tunneling rates in Fig. HJb). The SET current 
is higher for n = 1 than for n — [sec Fig.QJc)], inducing 



a higher excess heating power for the superconductor at 
n = l. However, even with the offset in the distributions 
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FIG. 4: (a) AsJ^t distributions obtained at T N = 142 mK with 
different drive frequencies, (b) Distributions of an individual 
measurement for forward and backward processes, (c) Bare 
entropy distributions for the measurements at Tat = 142 mK. 




FIG. 5: (a) Test of the detailed fluctuation relation, Eq. ([6|, 
for the measured trajectory entropy distributions. Despite 
the asymmetry of forward and backward processes due to de- 
tector back-action, the relation is satisfied, (b) The detailed 
fluctuation relation for bare entropy production of the forward 
processes. The solid lines in (a) and (b) show the expected 
dependence given by Eq. (JSJ). 
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in Fig.[4jb), they obey the detailed fluctuation theorem, 
as shown in Fig. [5] 

To summarize, we have extracted the distributions of 
work, bare entropy production, and coarse-grained en- 
tropy production for a sinusoidal drive protocol in a 
single-electron box. Due to the thermal non-equilibrium 
caused by the overheating of the S island, the work and 
entropy distributions no longer coincide. As a conse- 
quence, the Jarzynski equality is no longer applicable, 
but the integral and detailed fluctuation relations for en- 
tropy production are shown to be valid. We also find 
that the two different measures of entropy arc related by 
a universal formula. 
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Maisi, Frank Hckking, and Simonc Gasparinetti for use- 
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SUPPLEMENTARY MATERIAL Since the environment function satisfies detailed balance, 

P{E E )/ 'P(-Ee) = e^ BEB , and the fermi function satis- 
fies eP*"B E »/sf N/ s(E N/s ) = f N /s(-E N/s ), one obtains 
7±{E S ,Q,E E ) = -f T (E s ,Q,E E ), and with Eq. 



RELATION BETWEEN BARE AND COARSE 
GRAINED ENTROPY 

In the following, we derive the equality 

(e" As ">) = e~ As ™ (7) 

for trajectories with a single tunneling event in a NIS 
single electron box. Let n be the net number of electrons 
tunneled from S to TV. An electron may either tunnel 
from the S island to the N island, or tunnel from N-lead 
to the S-lead, adding or substracting one to the param- 
eter n. Let the tunneling events of the former type be 
denoted with '+', and of the latter with '— '. We consider 
a box at an electromagnetic environment at temperature 
T E . Upon a tunneling event, the environment absorbs or 
emits a photon with energy E E : the energy of the elec- 
tron in the superconducting lead Eg and the energy in 
the normal lead P/v then satisfy 



-As 



Eg — E 



N 



±(Q + E E ), 



(8) 



where Q is the dissipated heat upon the tunneling process 
0—^1, determined solely by the external control param- 
eter n g . The dimensionless medium entropy change of 
such an event is 



A4 = TPsEs ± PnE n ± p E E E , 



(9) 



where /3, = l/k E Ti denotes the inverse temperature. 

According to the theory for sequential tunneling in an 
electromagnetic environment Q, the tunneling rates for 
transitions Tq^i = T + and Ti^o = T- arc 



r±(Q) =J dE sJ dE El± (E s ,Q,E E ); 



-Y±{E S ,Q,E E ) 



e 2 R T 



N s (E s )f s (±Eg)x 



(10) 



P(±E E )f N ( T E N ), 



where Rt is the tunneling resistance, Ng{E) = 
is the normalized BCS super 



Re(j£|/V(£ 2 -A2)j 
conductor density of states with a superconductor en- 
ergy gap A, P{E E ) is the probability for the environ- 
ment to absorb the energy Q E , and /jv/sC-Ejv/s) = 
(1 + cxp(p N /gE N / s ))~ 1 is the fermi function of the 
N/S lead, giving the probability for an electron to oc- 
cupy the energy level E N /g. The conditional probability 
for the energy parameters to be exactly Eg and E E is 
P{E Sl E E | Q,n^n±l)= 1± (E s ,Q,E E )/T ± (Q). 
Left hand side of Eq. becomes 



-As 



} 



dEs 



TOl) the average is then 



r T (Q) 
r±(Q) 



-As 



(12) 



EXTRACTION OF TUNNELING RATES FROM 
STATE PROBABILITIES 

The tunneling rates Ti^o and ro_>i are obtained from 
the master equation for a two state system. At any given 
time instant t, the system has a probability Pi to occupy 
the charge state n = 1. As the charge state must be 
either n — or n — 1, the occupation probability for 
n = is Pq = 1 — Pi . The master equation is then 



Pi 



-IW,K(t))Pi + IW(n s (i))(l - Pi)- (13) 



In order to solve the tunneling rates as a function of n g , 
the occupation probability is calculated for both forward 
n~^(t) and reverse drives n^~(t). These satisfy n^(t) = 
rig~* (r — t), and by a change of variable t' = T—t, (t) = 
n^(t'), two equations arc obtained: 

A" 4 = -iW<(t))pr + r _>i(n^(t))(i - /T); 

-pt = -iw<(t'))pr +iwv(t'))(i - pt)- 

(14) 

The rates are then solved as 

tM 1 -p->) +p->-(i -p<-) 



ri_rt(n^(*)) 

iW(n^)) 



p<- -p- 

P^P^- + P<-P->- 



(15) 



P< 



P- 



and p<_ are obtained from the measurements for a 
time interval t...t + At by averaging n over the ensemble 
of process repetitions. The obtained distributions are in- 
serted in Eq. (fT5|) to obtain the rates. As shown in Fig. 2 
(b), Eq. (fTU|) describes the extracted tunneling rates well. 
The direct effect of the environment is negligible, and we 
take the limit of weak environment, P(E E ) = 8(E E ). 
By fitting Eq. (fT0|) to the measured rates, we obtain 
R T ~ 1.7 Mfi, A ~ 224 ^eV, and E c ^ 162 ^cV. Here, 
Tjv is assumed to be the temperature of the cryostat, 
while Tg is obtained for each measurement separately as 
listed in Table 1. 



FABRICATION METHODS 

The sample was fabricated by the standard shadow 
evaporation technique 0). The superconducting struc- 
tures are aluminium with a thickness of ~ 25 nm. The 
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tunnel barriers are formed by exposing the aluminium 
to oxygen, oxidizing its surface into an insulating alu- 
minium oxide layer. The normal metal is copper with a 
thickness of ~ 30 nm. 
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